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Abstract 

In view of the current discussion on the subject, an effort is made to 
show very accurately both analytically and numerically how the Drude 
dispersion model gives consistent results for the Casimir free energy 
at low temperatures. Specifically, for the free energy near T = we 
find the leading term proportional to T 2 and the next-to-leading term 
proportional to T 5 / 2 . These terms give rise to zero Casimir entropy 
as T — > 0, and is thus in accordance with Nernst's theorem. 
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1 Introduction 



In recent years there has been a lively discussion about the thermodynamic 
consistency of the expression for the Casimir pressure at finite temperature T. 
The problem gets accentuated at low values of T, where one has to satisfy the 
Nernst theorem saying that S = —dF/dT goes to zero as T — > 0. (Here S is 
the entropy and F the free energy, both referring to unit plate area.) What we 
shall consider in the following is the standard Casimir configuration, implying 
two semi-infinite homogeneous metallic media separated by a vacuum gap of 
width a. We take the two media to be identical, and assume that they are 
nonmagnetic with a frequency dependent relative permittivity e(u). Spatial 
dispersion is neglected. The two surfaces lying at z = and z = a are 
assumed to be perfectly plane, and to be of infinite extent. 

A central ingredient in the discussion of the thermodynamic consistency 
of calculated results for the Casimir attractive force between real materials 
is the form of dispersion relation used as input in the conventional Lifshitz 
formula. A very useful dispersion relation - the one that in our opinion is 
by far the most preferable one amongst simple dispersion relations for real 
systems at arbitrary frequencies - is the Drude expression 



Here us = i(, u p is the plasma frequency, and v is the relaxation frequency (we 
use the same notation as in Ref. [I]). The plasma wavelength is \ p = 2ttc/u p . 
For gold, the substance that we shall focus on in the following, we use 



In Ref. [T], we employed the values u p = 9.0 eV, v = 35 meV, what amounts 
roughly to a difference on the 1% level. The exact determination of Drude 
parameters is a non-trivial matter as discussed in [2]. Using a slightly dif- 
ferent set of Drude parameters will shift our numerical results slightly, but 
does not alter any of our conclusions. 

When comparing with experimental values it turns out that the Drude 
relation fits optical data very accurately for ( < 2 x 10 15 rad/s [31 Hj. In this 
connection we should bear in mind the following fact (cf. also the discussion 
in Ref. [5]): There exist no measurements of the permittivity at very low 
frequencies. What is available, is a series of measurements of the imaginary 




(1) 



u p = 9.03 eV, v = 34.5 meV, X p = 137.4 nm. 



(2) 
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part e"(u) of the complex permittivity e(u) = s'{u) +ie"(u). The Kramers- 
Kronig relations then permit us to calculate the real part e'(uj), and thus 
the complete e(u) is known. Permittivity data kindly supplied by Astrid 
Lambrecht cover a very large frequency region, from 1.5 x 10 11 rad/s to 
1.5 x 10 18 rad/s. From these data, the relaxation frequency v can be derived. 
As mentioned, from comparison with experimental data it turns out that v 
can be taken to be constant to a good accuracy, up to about ( = 2 x 10 15 
rad/s. For low frequencies the Drude relation yields the proper extrapolation 
down to ( = 0. 

A word of caution is called for, as regards the circumstance that permit- 
tivity measurements are done at room temperature in practice. The frequency 
v is in principle temperature dependent, and we do not know the value of 
v{T = 0) very accurately. It might seem natural here to invoke the Bloch- 
Griineisen formula for the temperature dependence of the electrical resistivity 
(cf. [6 J or also the discussion in Ref. pQ). From this, a form for v = v(T) 
can in principle be found. According to the formula, the value of v should 
go to zero as T — > 0. However, in practice this is not true. There are always 
impurities present, which make the value of v finite at T = [7]. The Bloch- 
Gruneisen formula, thus, is not followed in this limit. Mathematically, the 
important point is that 

C 2 [£«)-1]^0 as C^O. (3) 

This relation ensures that the zero-frequency transverse electric (TE) mode 
does not contribute to the Casimir force at all, as discussed in detail in 
Ref. [Tj . Strictly, the Drude parameters of Eq. fl2]) are valid at room-temperature, 
and v will take significantly smaller values for low temperatures. This affects 
our numerical results quantitatively, but not qualitatively; as long as v is 
nonzero, the TE part of the free energy vanishes at zero frequency which is 
the central point. 

The recent series of works on the Casimir effect by the present group 
of authors PQ El El HDl El [HI H21 US] - built upon the Lifshitz formula and 
the measured values of e(u) in combination with the Drude relation - have 
nowhere been found to run into conflict with basic thermodynamic principles. 
And there are other papers in agreement with ours: for instance, Jancovici 
and Samaj [33], and Buenzli and Martin [T3] considered the Casimir force 
between two plates in the high-temperature limit. They found that the linear 
dependence in T for the Casimir force to be reduced by a factor of two from 
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the behavior of of an ideal metal, this being a signal of the vanishing influence 
from the zero-frequency TE mode. (The first observation of the vanishing 
influence from this particular mode was made by Bostrom and Sernelius |16j.) 
Further support is found in the paper of Sernelius [17], who calculates the 
Casimir force taking spatial dispersion into account as well. It is found that 
at high temperatures and/or at large separations the force is reduced by the 
same factor of two compared with the ideal- metal result. 

There is no universal agreement on these issues, however. In a series of 
recent papers - cf., for instance, Refs. [T8J [191 I2Q1 [211 122] - it is argued that 
the Drude dispersion relation runs into trouble explaining recent experiments, 
and moreover comes into conflict with the Nernst theorem. These authors 
favor, instead of the Drude relation, the plasma relation 

e(<C) = 1 + ^f, (4) 

which corresponds to setting v = in Eq. (TjQ). (It should be noted that the 
expression (jlj) does not satisfy the condition (|3J).) 

An argument of the latter references is that omission of a zero frequency 
TE-mode would add a term linear in T to the free energy. This would violate 
the Nernst theorem as it would give a non-zero contribution to the entropy 
at T = 0. However, this argument is based upon use of the IM (idealized 
metal) model where e = oo for all (, or use of the plasma model §Qj where 
no relaxation is present. With v > the term in question is still linear 
away from T = 0, but the precise behavior as T —>■ has been less obvious. 
As argued in Ref. [1] the straight line should bend to become horizontal at 
T = 0. This was not verified in utmost detail, however; the numerical results 
of Ref. [TJ did not go sufficiently close to T = to show the behavior very 
distinctly, and the previous discussion and disagreement about violation of 
the Nernst theorem has accordingly continued. The main purpose of the 
present work is to investigate the issue more closely: we will show in detail, 
both analytically and numerically, how the Casimir energy behaves close to 
T = and by that show how it is consistent with the Nernst theorem. 

We shall not go into a study of experimental aspects in this paper. 
Rather, the objections referred to above makes it mandatory to reconsider the 
thermodynamics associated with the Drude relation anew, which brings us to 
the central theme of this paper. We will aim at showing, via a combination 
of analytical and numerical methods, how the basic theory sketched above 
(essentially the Drude theory) satisfies the Nernst theorem to a very high 



4 



accuracy. We consider this point to be important; a simple physical model 
of course cannot be permitted to run into conflict with thermodynamics. 

In the next section we show analytically, by using the Euler-Maclaurin 
formula, that the dominant contribution to the free Casimir energy F near 
T = is proportional to T 2 . We evaluate both this term and the leading 
correction term, which is proportional to T 5//2 . This implies that the entropy 
S = —dFj dT tends to zero as T — > 0, in accordance with the Nernst theorem. 
In Sect. 3 we calculate F numerically, and find agreement with the previous 
analytical result to a very high degree of accuracy. The results are illustrated 
in various figures. Thus we can conclude that the Drude ansatz does not run 
into conflict with thermodynamics at all. 

Readers interested in recent reviews on the Casimir effect may consult 
Milton's bookEI [23], and several review articles [2H [2S1 [26J, [27]. A great deal 
of recent information can also be found in the special issues of J. Phys. A: 
Math. Gen. [28] and of New J. Phys. [29]. 

As mentioned above, we shall not be concerned with comparison between 
theory and experiment in the present paper. We mention, though, the recent 
experiment of Obrecht et al. |30j, which seems to report the first accurate 
measurement of thermal Casimir effects. The experiment is important, but 
it lies outside the scope of the present investigation since it deals with non- 
uniformly heated systems. 

Finally, we mention the special variant of the thermal Casimir problem 
consisting in studying, instead of a metal, a semiconductor endowed with a 
small but finite conductivity at zero frequency [31, 32J. According to the 
authors of these references this situation implies an interchange of the roles 
of the TE and TM (transverse magnetic) modes, as compared with the case 
of a metal. Namely, within an idealized approach, they find that the TM 
reflection coefficient gets a discontinuity at ( = 0, implying in turn an ap- 
parent conflict with the Nernst theorem. The problem is interesting, and we 
hope to return to it in a later paper. 

5 It may be mentioned for completeness that this book from 2001 was written from the 
standpoint of the "classical" IM model. 
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2 Analytical approach: Casimir free energy 
near T = for real metals 



As mentioned, we use the same notation as in Ref. pQ. We will evaluate the 
leading T-dependence of the Casimir free energy near T = for metals using 
the Drude relation ([T]), 



D 



with D = ujp/u. The free energy is given by expression (3.4a) in pQ as 

PF = ^ E ' /Kl - A™) + ln(l - A TE )]gdg (6) 
Ztx — „ J 

m=0 c/c 

where /3 = 1/A;T and 

A™ = Ae- 2(?a , 
A TE = Be- 2qa . 

The prime on the summation sign means that the case m = is taken with 
half weight. The coefficients A and B are the squared Fresnel reflection 
coefficients for the two media and are given by 

A=(^) 2 , B=( S -^)\ (7) 
\s + ep J \s+pj 

with 

s = ^Je-l+p 2 , p = — . 

Here a is the plate separation, e the relative permittivity, c the velocity of 
light in vacuum, and ( is the Matsubara frequency given by 

, 2ttA; _ 

C = —mT. (8) 

Note that the quantities of A, B, s,p and ( all depend upon the summation 
variable m. (Units c = h — k — 1 are not used.) The term of interest is the 
TE mode, since this is the term that gives rise to the controversy about the 
Nernst theorem. 
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With the small £ dependence of Eq. (J5]) the .B has a scaling form such 
that it can be expressed in terms of one variable. So by introducing a new 
variable x to replace q the (^-dependence can be removed fully 

With this we have 



'VT+ 



B=[ V : ' X 2 * ) =(v / TT^-x) 4 , (10) 
and the TE free energy expression can be written as 

oo ' 

f3F TE = Cj2 9(m), (11) 



m=0 

with 

g(m) = m 

where 



roo 

/ aTn(l- Be- y )dx, (12) 



The small T-behavior can now be obtained by use of the Euler-Maclaurin 
summation formula 

f>(n) = ]g{u) du - l</(0) + ±g"'(0) - ... (14) 
One easily sees that 

poo 

g'(0) = / xln (I -B)dx. (15) 
Jo 

This integral can be performed analytically. First introduce a new variable 
x = sinh t with dx = cosh t dt 

oo oo 

g '(0) = J sinht cosht ln(l - e~ u ) dt = ^J {e 2t - e~ 2t ) ln(l - e" 4t ) dt. (16) 
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Next we substitute e 21 = u for which — 2e 21 dt = du and use partial inte- 
gration to obtain 



9(0) = ~J (l-^)ln(l-n 2 )^ 
o 

= -i + it) 2 ln(l + m) + (1 - u) 2 ln(l - u)] - 2u) 

o U Jo 

= --(21n2- 1) « -0.09657. (17) 

At T = the free energy is determined by the integral in Eq. ( il4"|) (besides 
contribution from the TM-mode). For small T — * the deviation from the 
T = value is thus {(3 = 1/kT) 



-->) 



C 

'21n2-l). (18) 



48/? 



This result was presented by Milton at the QFEXT03 Workshop, and is given 
as Eq. (22) in Ref. [8], or Eq. (4.9) in Ref. p]. 

It can be noted that AF TE is independent of the plate separation a and 
can thus be valid only sufficiently close to T = for a given a such that 
AF TE <C F TE ~ he/ a 3 . Evaluation of the next term as given by result (|30|) 
below verifies this. The leading term of AF TE dominates its next term only 
when the dimensionless quantity aC 1 ! 2 < 1 or a 2 kT hue 2 /uj 2 . A conse- 
quence of this is that for increasing a the temperature interval where ( Tl8l) is 
valid decreases rapidly, and AF TE becomes more and more negligible com- 
pared to F TE since then AF TE ~ C//3 = u 2 (kT) 2 /(c 2 hu) < c 2 hu/(uj 2 a 4 ). 
(Thus in the present case AF TE F TE provided a > cv/uj 2 ~ 0.1 nm with 
a 2 kT << hvc 2 /uj 2 fulfilled.) In the limit a — > oo the IM (e = oo), but now 
with the linear term included, is recovered. But the latter does not violate 
the Nernst theorem as long as a is finite, and for a = oo there is no Casimir 
free energy anywaj{^]. 

Inserting the value for g'(0) and the values huj p = 9.03 eV and hv = 
34.5 meV for Au we find with C given by (IIBj) (h = 1.0545 • 10" 34 Js, k = 



6 It may be noted that for a/(Tic(3) <C 1 an a-independent contribution to the free 
energy was found also in the "classical" Casimir theory for metals at low temperature. 



See Eq. (3.38) in Ref. [23] or Eq. (3.24) in Ref. [J. 
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1.381 • 1(T 23 J/K, c = 2.998 • 10 8 m/s) 

AF TE = ^2 with ^ ^ 5>gl . 1Q -13 J/( m 2 K 2 ). (19) 



It turns out that Eq. ([TBI holds only very close to T = (i. e. T 
0.01 K), but there will be a leading correction that we can obtain with good 
accuracy. Expanding g(m) in powers of m one notes that half integer powers 
will occur. Thus formula (|14j) is not quite valid as g"'{0) and higher order 
derivatives will diverge. However, this problem can be avoided since the 
formula can be applied to summation starting at m — p (p — 1, 2, 3, • ■ ■). 
Thus we have 

p-i 



oo p-± ~ 111 

E #( n ) = E + y 9(u) du + -<?(p) - — g\p) + ^"'(p) + ■ ■ • 

/ du + S- -g\p) + — p"'(p) + ■ ■ • (20) 



n=0 n=0 

oo 



where 



p-i 1 p 

n=0 o 



For a power term g a (m) = m" (for small m) we have 

p -. 

£(p) = <nr-\ flip) = a(a - l)(<r - 2)p°~\ f g a {u) du = p°+\(22) 

o 

With the choice p = 1 we get 

S = SJl) = - — . (23) 

One may note that Si(l) = as should be expected. The power of key 
interest here will be a = 3/2 by which S a — 1/10 , and thus 

AS a (l) = S a {l) - -^(1) + ^5^(1) ~ -0.02552 « -0.204^(1). (24) 

When other terms are neglected the error can be estimated from the next 
term in the series (provided the sum of higher order terms converge) 
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4.65- 10" 



(25) 



which is only about 0.2% of the value (l24"j) . For increasing values of p 
the error goes further down rapidly. (Instead of result (1241) p = 2 yields 
~ —0.02549.) Thus ff24|) with p = 1 is a good estimate for the sum of 
interest . 

With g(m) given by (fT2"j) we can expand to leading order in m or y (oc 
T 1 / 2 ). To this leading order the limit of integration xq (oc T 1 / 2 ) can be put 
equal to zero since the integrand vanishes for x = 0. We find 



g{m) 



mg'(0) + myl + 



Now from (ED and (EE 



y = 2av / 27Ca;m 1/2 , 
with C given by (fl3|) . Thus the derivative becomes 

g\ m ) = g '(Q) + 3aV2^CIm 1/2 + ■ 



(26) 
(27) 



(28) 



where g'(0) is given by (fTTj) and (with substitution a; = sinht as in ffTB"]) ) 



" o; 2 B , 1 " 
& = ^ 



-t 



-3i\ 



dt 



o 







1 

12' 



(29) 



Now we will use Eq. (121]) for the most simple case m = p = 1 for which 
S a (l) is given by (]2"3"]) . and with exponent a equal to 1 and 3/2 for the 
two separate terms in (I26p the corresponding values of S a (l) are and 1/10 
respectively. Thus with flM} and ([28]) we have 12 Aft (1) = -#1(0) = -#'(0) 
and 12 AS 3/2 (1) = -0.204 g' 3/2 (l) = -0.204 • 3aV2^CI. So with JTBJ the 
change in free energy becomes 



AF 



TE 



2[A5 1 (1) + Aft/ 2 (1) + 



(3 



-—g'{0) 



0.204-^^ + 
-12^(0) 



(30) 
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which gives the small T correction to the result (TT8|) . Again inserting the 
previous values for u p and v we obtain with plate separation a = 1000 nm 

AF TE = dT 2 [l -C 2 T^ 2 + ■■■] with C 2 « 3.03 K- 1 / 2 . (31) 

This result can only be valid for very small T as it otherwise becomes 
negative already for T slightly larger than 0.1K. To avoid this one may 
instead write 

AF * = T%vr < 32 > 

as the theoretical or analytical result for AF TE for small T. This Pade ap- 
proximant form (which is equivalent to Eq. (l3~Tj) with respect to the first 
two terms), turns out to be convenient for comparison with numerical eval- 
uations. If the corresponding numerical result for AF TE is AFj u ^ one can 
evaluate the ratio 

AF TE — AF TE 

rj th num fQQ^ 



and consider the limit T — > for which the limiting value should be R = 
(cf. more details in Appendix B). 



3 Numerical calculation of free energy at low 
temperatures 

We have calculated the free energy F TE as a function of temperature given by 
([6]) for two gold half-spaces separated by a vacuum gap of width a = 1.0 /zm. 
This would be a typical experimental situation where the influence from 
the finite temperature is large, about 15% increase in the magnitude of the 
Casimir free energy PQ, and a corresponding decrease of the Casimir force 
according to our theory with no contribution from the TE mode at zero 
frequency. The calculations use the permittivity data for gold, received from 
Astrid Lambrecht as mentioned. At T = the free energy is calculated 
numerically as a double integral rather than a sum of integrals using a two- 
dimensional version of Simpson's method. 

The vanishing of the zero frequency mode is intimately connected with 
the behavior of the reflection coefficient B at vanishing frequency. According 
to the Drude model (or any model satisfying (j3J)) the TE mode reflection 
coefficient tends to zero as £ —>■ 0. To illuminate this point, we have plotted 
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Figure 1: Squared reflection coefficients A and B of the metal interfaces 
as defined in ([7]) for the TM and TE modes, as a function of £/c and the 
transverse momentum k±. a) A for the TM mode, b) B for the TE mode, c) 
B for kj_ and £ close to zero. 



12 



Plate material : Gold 




Temperature (K) 



Figure 2: Numerical evaluation of the free energy ([6]) between two gold half- 
spaces as a function of temperature. The inset gives details for low T. 



A and B as a function of imaginary frequency and transverse momentum, k±, 
for an interface between gold and vacuum in Fig. [TJ In part c) of this figure, 
we clearly see how B vanishes smoothly when Q — > for k± ^ consistent 
with Maxwell's equations of electrodynamics. However, the coefficient A in 
Fig. UK) for the TM mode equals 1 for all k± as ( — > 0, as is also the case for 
for all ( for an ideal metal. In the latter limit we also have B — > 1, but for 
( = 0, B remains zero. 

For ideal or non-ideal metals it is well known that the temperature cor- 
rection for the TM mode behaves as T 4 . Thus it does not add to the T 2 and 
rp5/2 terms that we find from the TE mode . 

By direct numerical integration and lengthy summations independent of 
the analytic derivations made in Section [2J we obtain the free energy numer- 
ically. Figure [2] shows the free energy versus temperature up to 800 K, while 
the inset shows details of the parabolic shape close to T = 0. First of all 
the figure shows the controversial decrease of the magnitude of the Casimir 
free energy and thus also the related decrease of the Casimir force up to a 
certain temperature. Secondly, the inset shows that the tangent of the curve 
is horizontal at T = as predicted. This implies that the entropy due to 
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Figure 3: Plot of the ratio R defined in Eq. ( 1331) . 



the Casimir force is indeed zero at T = 0. Thus the Nernst theorem is not 
violated when using the realistic Drude dispersion model. In contrast, it 
would be violated if a TE-term were added for £ = 0. This conclusion, as 
mentioned above, is clearly in contrast to that presented in various earlier 
works [TH [HI [201 ETJ [22] . The reader should note that the dependency of v 
on temperature has been neglected in figure [21 

Now there are always some uncertainties connected with numerical calcu- 
lations. Also the analytic derivation has some uncertainties , e.g., concerning 
proper use of the Euler-Maclaurin formula, and concerning convergence and 
neglect of of higher order terms. In Figj3] we have therefore made a more 
accurate and much more sensitive test of the behavior near T = comparing 
the analytic result with the numerical one, by plotting the ratio R defined 
in Eq. (J33l) . We see that R when extrapolated approaches zero linearly with 
a finite slope (when taking the curvature of the plot into account). Thus, 
with high accuracy we find full agreement concerning the T 2 and T 5 / 2 terms 
in the free energy and their coefficients C\ and C2. As the number of terms 
to be summed numerically increases rapidly when T = is approached, our 
evaluations were terminated at T = 0.008 K. The extrapolated value R = 
for T = means that the coefficient C\ is correct while the finite slope means 
that C% is correct too within numerical uncertainties. Also if a more dom- 
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inant power were present, R would diverge at T = 0. The finite slope of 
R at T = means that the next term in AF TE is of higher order (see also 
Appendix B for details). 
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A Alternative derivation by expansion of g(m) 

As a variant of the analytic approach, let us show how the essential depen- 
dence of the free energy F TE on T near T = also can be recovered by 
making use of complex integration. We begin with the TE expression 

PF TE = CJ2™ xln(l - Be~ ax )dx, (Al) 

m=0 Jx ° 

where 

J 1 

C = -^7, a = 2aV2nCm. (A2) 
pave 1 

The essential step now is to expand the logarithm to the first term, 

°° /"OO 

(3F TE = -C Y m / xBe~ ax dx, (A3) 
Jo 



m=l 



(m = does not contribute to the sum (1A1H ) . where we have also replaced 
the lower limit xq by zero. We next use the formula [31] 

1 rc+ioo 

e~ ax = — ;/ ds(ax)- s T(s), c> 0, (A4) 
where T(s) is the gamma function. The summation over m is easily done, 

oo 

E ™ W2 = C(f - 1)- (A5) 

m=l Z 

Here ( is the Riemann zeta function, defined by the analytic continuation of 

C(*) = f7T/ -r-rdt (A6) 

1 (s) Jo e — 1 
15 



for Re(s) < 1 [33J. As is known, ((s) is one- valued everywhere except for 
s = 1, where it has a single pole with residue 1. As T(s) has simple poles 
at s — —n with residue (— l) n /n\, n — 0, 1, 2, ... we get, by taking c > 4 and 
closing the contour of integral flA4j) as a large semicircle on the left, 



TE 



-C 



Bx 



r(4) 



(2ax) 4 (2ttC) : 



((-1) -2axV2^C ((--) + ... dx. 

(A7) 

The first term in (1A7|) diverges, the reason being that we have replaced the 
lower limit xq by zero. This term is independent of T, and can thus be omitted 
in the present context. Thus using ((-1) = -1/12, C(~3/2) = -0.025485, 
together with the integrals J °° Bxdx = 1/12 and / °° Bx 2 dx = 8/105 gives 
the temperature-dependent free energy to first order in B. 

However, this is easily extended to arbitrary order in B by expansion 
of the integrand in flAlj) since (B e~ ax ) n = B n {l — anx + |(cmx) 2 + ■■•). 
Thus by summation xB — > xY^=\B n /n = — xln(l — B) and x(Bx) — > 
xY^=i{nB n / n)x = x 2 B/(l — B). By this the above integrals become those 
of Eqs. ffTBT) and ([2"5|) . and for the temperature dependent part of the free 
energy we thus get 



AF 



TE 



C 
1 



c 

1 



V(0)C(-l)-2av / 27C/Cl --) + 
3\ 3aV2^C 



12 



i 



'C 



-12g'(0) 



+ 



(A8) 



fully consistent with result fl30|) . Note that £(— 3/2) is close to the approx- 
imate value ff24l) and even closer to the more accurate value ff25|) . Thus 
we have reason to consider £(—3/2) to be the exact result for the Euler- 
Maclaurin expansion performed in Sec. [2j 



B Remarks on the quantity R 

Let us explain in some more detail the interpretation of the quantity R, 
defined in Eq. fl33l) . This quantity gives the relative difference between the 
temperature dependent theoretical free energy AF^ h E , having the form fl32|) . 
and the temperature dependent numerical free energy AF^ E m calculated from 
Eq. ([6]). (As mentioned in Sec. 3, the TM mode behaves as T 4 and is thus 
negligible near T = 0.) 
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Let us assume that has the same form (1321) as before, with coeffi- 

cients C\ and C 2 , and that /S.F^ m has the form 

^IL = A(T 2 - D 2 T 5 / 2 + D 3 T 3 + ■••), (Bl) 

with calculated values for the coefficients D%, D 2 , and D%. Then, 

(B2) 

If C\ = Di and C 2 = D 2 , we see that R is zero at T = and is linear in 
T for low T. From Fig. 3 we see that the fit is perfect insofar as it may be 
determined from the graph. A constant term would have caused a nonzero 
value at T = 0, and a nonzero T 1//2 term would have caused a vertical slope 
near T = 0. None of these effects are perceivable within the numerical 
accuracy, from which we must conclude that C\ and C 2 are correct within 
numerical accuracy. 
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